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Abstract 

We show that particles can spht only when their group velocity exceeds their phase velocity. 
I In this sense the splitting process is the quantum analog of the modulational instability 
I in anomalous dispersive media. In the case of a neutrino whose dispersion relation con- 
^ ! tains a subdominant Lorentz- violating correction of the form aP^, the neutrino will decay 
into two neutrinos and an antineutrino at a rate proportional to a^G^E'^'^^'' . Unlike the 
Cohen-Glashow instability, the splitting instability exists even if all particles involved in the 
^ interaction have the same dispersion relations at the relevant energy scales. We show that 
I this instability leads to strong constraints even if the energy ii^ is a function of both the mo- 
mentum P and also of the background density p, for example we show that it alone would 
have been sufficient to eliminate any model of the MINOS/OPERA velocity anomaly which 
modifies the neutrino dispersion relation while leaving those of other particles intact. 
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1 Motivation 



The velocity anomaly reported by OPERA [T] was the result of experimental error [2] . 
The frenzied activity following OPERA's announcement is often criticized for its inconsis- 
tent attempts at profound advances in fundamental physics. However it also led to several 
noteable advances of general interest, in particular with regard to constraints on unusual 
dispersion relations. The most celebrated among these is that of Cohen and Glashow [3] 
who demonstrated that if at high energies the velocities of two species of particle asymptote 
to different values, then the fast particle will lose energy as it travels, transferring it into a 
succession of slow particles. Similarly it was demonstrated that the phase space for a slow 
particle to decay into fast particles is extremely limited [1]. 

These two new constraints were powerful enough to eliminate all of the models which 
had been formulated at that time, even those with exotic dispersion relations such as those 
of Refs. [SI E]- For example, consider the model of Ref. [S] in which the neutrino velocity is 
a power of the momentum P, with k > 0. As we will explain, when k > 0.5 these models 
were ruled experimentally by OPERA's short burst run [1] in which the energy varied by 
at least a factor of 5 but the velocity varied by less than a factor of 2. On the other hand 
k > 2.5 is required for consistency with supernova bounds [3 |8]. These models can also be 
ruled out theoretically, as is easily seen by extending Cohen and Glashow's calculation to 
the proposed dispersion relation. Such generalizations were considered in Refs. [HI [TOl [TT] , 
where it was seen that they only change the decay rate by a /c-dependent geometrical factor 
of order 1, not by nearly enough to avoid the Cohen-Glashow electron-positron nucleation 
process. As a result the neutrino energy decays appreciably before arriving at OPERA, in 
contradiction with the fact that the observed spectrum agreed well with simulations up to 
100 GeV [12] . 

The model of Ref [6] is consistent with the results of the short burst run, as the velocity 
is constant throughout the energy range probed. However the dispersion relation is nearly 
identical to that of Cohen and Glashow, with the only difference arising in a small portion of 
phase space at small momenta. Thus again an application of Cohen and Glashow's argument 
implies that in this model the neutrinos should have lost most of their energy before arriving 
at the OPERA detector, in contradiction with observations. 

Refs. [31 H] remain relevant even now that the velocity anomaly has disappeared because 
they provide general constraints on dispersion relations that differ from those of special 
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relativity. As is clearly explained in the first two paragraphs of Ref. [13], the dispersion 
relation of special relativity only applies to particles in the vacuum. In particular it applies 
to particles in a Lorentz-invariant vacuum. It therefore does not apply to many systems of 
interest, such as neutrinos traveling through the Earth or particles in the primordial plasma. 
If the acceleration of the universe is caused by an interacting dark energy field, then the 
expectation value of this field is not Lorentz-invariant, it depends on a universal time. As 
a result, at some scale the dispersion relations of all particles will deviate from the special 
relativistic form. Even if dark energy is simply caused by a cosmological constant, the FRLW 
solution of our universe is not Lorentz-invariant, and so one may expect a combination of 
Hubble and Planck scale corrections to the relativistic dispersion relations. 

More immediately, neutrinos traveling inside of media still exhibit a number of anomalies. 
For example there is a solar neutrino deficit at low energies [H] and both AMANDA and 
IceCube observe a surplus of horizontal neutrinos in ice [15]. LSND, whose baseline is 
dominated by Earth and iron shielding has seen an anomalous neutrino deficit [16] while 
KARMEN, whose baseline is dominated by air, has not [17]. Needless to say, these media 
are not Lorentz-invariant and so the neutrino dispersion relations inside of them, even within 
the context of the standard model [18], are not Lorentz-invariant. While it is not yet clear 
whether these anomahes are due to new physics, they motivate a better understanding of 
constraints on neutrino dispersion relations in general and in particular of neutrino dispersion 
relations in media. 

While the constraints [31 II] are quite powerful, eventually models of superluminal neu- 
trinos were created which avoided the assumptions implicit in these works and so avoided 
the constraints. The most obvious assumption is 4-momentum conservation, and discussions 
of escapes along those lines first appeared in Refs. [19] . It was also argued that one could 
escape by rendering other particles, such as the electron, superluminal [201 EI] after all the 
dispersion relations could well be density dependent [221 ESI 121] and all strong constraints on 
electron velocities appear to come from experiments in a vacuum [2l]. Another escape route 
is if, while traveling, the neutrinos are noninteracting because they are taking a short cut 
through another dimension [25], or converted into a sterile flavor [26]. While such models are 
no longer necessary to explain OPERA's anomaly, they may well resurface in attempts to 
explain the yet unexplained neutrino anomalies. This leads one to ask if further constraints 
are available which may eliminate some of these models. 
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2 Neutrino splitting 



2.1 What is neutrino splitting? 

In this note we will argue that neutrino splitting pTf 128] 

V — ^2z/ + z7 (2.1) 

provides a powerful constraint. Like the Cohen-Glashow and phase space constraints, the 
applications of this constraint extend far beyond OPERA. The Cohen-Glashow process is 
only allowed if two different kinds of particles have different dispersion relations. Therefore 
it is avoided in models with universal high energy dispersion relations such as plj 120] [21] . 
While qualitative features of such models were discussed in Ref. [20], a concrete model of this 
kind was proposed in Ref [2^, considering interface conditions between regions with different 
densities and demonstrating that fifth force constraints can be satisfied in those regions. The 
phase space constraint [21] uses the kinematics of the creation of the neutrinos to limit its 
propagation along the baseline. However these neutrinos are created in a vacuum whereas 
the baseline is almost entirely solid rock, therefore even a small amount of environmental 
dependence is sufficient to eliminate these constraints. In other words, the constraint of 
Refs. [1] are satisfied in any theory in which the neutrinos are subluminal in a vacuum, such 
as those of Refs [201 EI]. 

The neutrino splitting constraint is more robust. It depends only on the dispersion 
relations of the neutrinos themselves, so it cannot be avoided by modifying the properties of 
the other particles. Furthermore the splitting process can occur anywhere along the baseline, 
so unlike the phase space constraint of [1] which only constrains the dispersion relations at 
the point where the neutrinos are created, the splitting constraint constrains the dispersion 
relations everywhere that the neutrinos travel and so in the region in which the velocity 
anomaly is claimed. 

When can splitting occur? Consider a field 0. Whenever the low energy effective theory 
has an effective vertex of the form 0^ with k > 2 there is a process in which a quantum of 
the field (p can decay into k — 1 quanta of (p. If the field carries a conserved charge then 
charge conservation does not allow an effective vertex of the form (p^. However, a vertex 
which includes k copies of the field and k of the antifield is allowed. For example in the case 
of neutrinos the 4 neutrino Fermi vertex is allowed, which yields the desired process (12.11) . 



^If the conserved charge generates a cycUe group of finite order N then this remains true so long as N 
does not divide k. 
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2.2 Splitting kinematics 



When is splitting kinematically allowed? When there is a large phase space for splitting, 
particles do not propagate very far before losing all of their energy. To obtain the strongest 
possible bounds one is interested in the regime where splitting is near threshold. In this 
case a particle with an initial momentum P^°) in the x direction will decay into k particles 
with momenta P'-'^ which are nearly in the x direction. More concretely, if for simplicity we 
assume that the x direction is a principal direction of our dispersion relation, then the energy 
of the original particle and the energy E« of a given outgoing particle are approximately 
quadratic in the corresponding momenta 

= ^(Pf)), E^'^ = E{P^'^) + a(Pf (2.2) 

where a is a strictly positive constant, E{p) is a function and the y axis is chosen such that 
P*^'^ lies on the x — y plane. 

In this approximation, we will now demonstrate the following claim: 

Claim: Particle splitting is only kinematically allowed when the group velocity exceeds the 
phase velocity. More precisely, it is only allowed for particles with x-momentum Pi°^ such 
that, at some x-momentum p < 

em > m, ,2.3) 

op p 

First, we note that Eq. (12.21) implies that the energy of the iih particle is greater than 
or equal to E{Px^) with equality only when Py^"* = 0. This equality is only satisfied on a 
measure zero portion of the phase space and so such decays do not contribute to the decay 
amplitude. This reflects the familiar fact that if a particle can only decay into colinear 
particles then the instability is only marginal. Therefore, observable decays only occur when 
the energy E^"^^ of each particle is strictly greater than the corresponding function E{Px^) 

E^^ > E{P^^). (2.4) 

Now we are ready to prove the claim by contradiction. Assume that the group velocity 
is less than or equal to the phase velocity 

am < m (2.5) 

dp p 

for all momenta p < P^\ Dividing both sides by E{p) one obtains 



dlniEip)) 1 , , 

op p 
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Now one may calculate E^'-^ using the fundamental theorem of calculus and then Eq. (12.61) 

' (2.7) 



which can be exponentiated to yield the second inequality in 

The first inequality follows from Eqs. (12. 2p and (12. 4p . Now summing over all of the final 
particles the rightmost and leftmost expressions yield unity. Thus we have proven that 
1 > 1, a contradiction. 

This implies that the initial assumptions that there is finite phase space for a splitting 
process and that the group velocity never exceeds the phase velocity at p < Pi°^ are incom- 
patible. So splitting is only allowed at momenta greater than momenta at which the group 
velocity exceeds the phase velocity, at least when one is close enough to threshold that the 
nearly colinear approximation (12. 2 p may be applied. 

What about the converse, does a group velocity faster than the phase velocity imply that 
the splitting process is kinematically allowed? If the group velocity is faster than the phase 
velocity at some momentum then at momenta p infinitesimally below P^\ the energy 
E{p) will be less that p^^- Let Pm be the lowest such momentum, so that for Pm < p < P^^ 

E{p) < (2.9) 

Then splitting is kinematically allowed so long as all of the decay products have p > Pm- The 
problem is that if Pm > -Pi^V^ then conservation of momentum does not allow all neutrinos 
to have momenta p > Pm, since the total momentum must be Pj'^^ 

Therefore we learn that if (12. 9 p holds for all p such that Px^^ /3 < p < Px^^ then the 
splitting is kinematically allowed. However even if it does not then splitting reactions of the 
kind 

au — >{a + b)u + bV (2.10) 

are allowed where (a + 2b) Pm < clPx^\ For any choice of b, Pm and Px^\ a sufficiently large 
value of a always satisfies this constraint. This means that a beam of neutrinos can always 
split if the group velocity exceeds the phase velocity, but the splitting rate may be highly 
suppressed if the phase velocity is only exceeded for a small portion of the set of momenta 
beneath Px^\ 
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2.3 Splitting rate 



Below we will calculate the splitting rate precisely in a particular model. In this subsection 
we will provide a rough estimate of the splitting rate for neutrinos. The rate is quite similar to 
that of the decay rate calculated by Cohen and Glashow [3], as the relevant process is virtually 
identical, one need only replace the electron-positron pair with a neutrino-antineutrino pair. 
There are still two vertices, leading to a factor of Gp in the amplitude and so G'jp in the 
decay rate. 

As we consider more general dispersion relations, there is not necessarily any limiting 
velocity and so the phase space computation is in general more complicated. There are still 
3 final particles, which have 12 4-momentum components. The conservation of 4 momentum 
yields 4 constraints and the mass shell conditions on the 3 particles yield 3 more, leaving 
12-4-3=5 undetermined quantities. Thus the phase space for the decay is 5-dimensional. 
Two of these five variables describe the fractional distribution of momenta between the 3 
final particles and 3 describe the orthogonal momenta. The orthogonal momenta, as we 
have described, would be constrained to vanish if the dispersion relation were linear. Thus 
these 3 variables live inside of a region whose radius is proportional to the extra phase space 
provided by the nonlinearity of the dispersion relations. 

The size of this extra phase space depends upon the choice of neutrino dispersion relation. 
For example, suppose one is interested in splitting at high energies and if at these high 
energies the dispersion relation is effectively E = P + aP^ with k > 1, where a has energy 
dimension [1 — A;]. Furthermore consider energies low enough so that aP^ << P, so the 
relativistic term is still dominant. Now the extra energy available is proportional to aE'' 
and so the phase space contribution to the decay rate is proportional to a^. 

Putting these results together, the decay rate is proportional to G'jpa^ which has energy 
dimension [—1 — 3A;]. The decay rate should have units of energy, therefore by dimensional 
analysis the decay rate is proportional to 



In the remainder of this note we will illustrate the power of neutrino splitting constraints 
by using them to demonstrate that the neutrino velocity anomaly could not have been caused 
by any modified dispersion relation for the neutrino alone, in other words a model consistent 
with opera's claim would have to modify the dispersion relations of other particles, which 



r ~ gWe^^'''. 




3 Splitting and OPERA 
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would be very difficult to achieve consistently with experimental constraints. 

Following Refs. [23t 120] we will consider dispersion relations E{P,p) which are general 
functions of the norm P of the 3-momentum and background density p. An arbitrary de- 
pendence on the density p is included in order to avoid friction with bounds on neutrino 
superluminalities from SN1987A [3 |8], which have plagued previous attempts at the deter- 
mination of such dispersion relations [9l [TOl HI]. The SN1987A bound may be satisfied by 
chosing the p-dependence such that at small p neutrinos are subluminal. Alternately the 
density-independent case, which includes for example the models of Refs. O IH] is included in 
our analysis as the special case in which the p dependence is trivial. In what follows we will 
make no assumptions concerning the p-dependence, our arguments will be completely gen- 
eral. However we will assert that the dispersion relations of all particles except for neutrinos 
are as in the Standard Model. 

To reduce the number of dispersion relations that we must consider, we use two tricks. 
First of all, we restrict our attention to experiments inside of the Earth's crust. Here the 
density is more or less constant. In a region of constant density po, the density-dependence 
of the dispersion relation is irrelevant and so for these experiments one need only consider 
a function in one variable E{P) = E{P,pq). Second, we impose a kind of monotonicity 
condition on E{P). In practice we impose that it has at most one inflection point by 
expanding to second order in P within a given interval, however this assumption can be 
weakened to allow as many as 10 inflection points without changing our conclusions. If, on 
the contrary, there were numerous small windows in P at which the neutrino is extremely 
subluminal, but these windows were so small than no neutrino yet observed by OPERA lies 
inside of one, then our conclusions would be evaded. 

With these restrictions on the function E{P) we apply three experimental constraints. 
The ffist constraint arises from the energies of the neutrinos observed in the 2008 and 2009 
OPERA runs and reported in Ref. [12]. This spectrum has two features which will be relevant 
here. First of all, it has a very long high energy tail, with 4 percent of events above 100 GeV. 
Second of all, it ffi expectations from simulations remarkably, even at energies of up to 100 
GeV. The high energies are relevant because the effects which will be considered are much 
stronger at high energy, where more phase space is available, and so they imply stronger 
bounds. This constrains the energy loss of the neutrinos during propagation, which places 
constraints on neutrino splitting and therefore as we will see on the excess of the group 
velocity over the phase velocity. 

The second experimental constraint comes from OPERA's short extraction run, which 
was able to unambiguously determine, up to whatever systematic errors there may be, the 
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anticipated arrival times of 20 neutrinos [T]. All of them arrived earlier than they would 
have had they traveled at the speed of light. Their velocities showed fractional excesses 
e = (f — c)/c of between 1.7 and 3.9 x 10^^. The energies of these neutrinos are not 
known, indeed 14 interacted outside of the detector. Nonetheless, they can be determined 
statistically using the energy distribution reported in Ref . [12] . The velocities determine the 
derivative of the dispersion relation in the momentum regime probed by OPERA. 

The final constraint that we will use is from the ICARUS experiment, whose detector is 
subjected to the same neutrinos as its cavemate OPERA. It did not observe any high energy, 
nearly collinear, e~e"'" pairs [29]. The resulting bound can be strengthened and extended to 
higher energies by results from the NOMAD experiment [30] and ICECUBE , but there 
is not yet evidence for neutrino superluminality at these energy scales and so this will not 
strengthen our claim. Generalizing Cohen and Glashow's results [3] for this process to a 
more general dispersion relation, this places a very strong bound on E{P) — P. 

These three constraints on the neutrino dispersion relation at OPERA energiej^ force 
it, in this energy range, to respectively not be too concave E'{P) < E{P)/P (neutrino 
splitting), to have a slope of 1 + e (neutrino superluminality) and to approximately lie in the 
triangle E{P) < P (Cohen-Glashow Bremsstrahlung) . The purpose of this note is to make 
these bounds quantitative, and so to argue that no such function exists even if the bounds 
are chosen very conservatively. Therefore a model which explains OPERA's superluminality 
claim must do more than just modify the dispersion relation of the neutrino. In the conclusion 
we will briefly mention some such additional features that have been considered in the 
literature which may render such a model consistent. 

In Sec. Hlwe will discuss the three experimental constraints that we will need to derive this 
argument. To derive the most robust bound, we will only assume neutrino superluminality in 
the energy regime in which it may be inferred from OPERA data. In Sec.[5]we will determine 
this regime and so use the bounds to determine the features of a dispersion relation which 
are necessary to satisfy the three constraints. We will then consider the dispersion relation 
which minimizes the neutrino splitting probability, and show in Sec. [6] that nonetheless it 
yields so much splitting as to be incompatible with OPERA's results. 

^We do not have these bounds at lower energies, but we do not need them. We consider only neutrino 
splitting in which all of the resulting neutrinos are within this energy range. The inclusion of decays with 
products below this energy range could only increase the decay rate, and so strengthen our result. 
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4 The Three Constraints 



There are many functions of the form E{P,p). We will greatly simplify the analysis by 
using only constraints at a fixed value of p. More precisely, we will consider only experiments 
in the Earth's crust, approximating its density to be constant. In particular, as most of the 
travel time of neutrinos from SN1987A was in much lower density media, we will not use 
these to fix any constraint. Fixing p, we arrive at a function only of momentum Ep{P). 
As we will be interested in experiments in which the momenta are roughly parallel to the 
surface of the Earth, we will ignore the potential directional dependence and only consider 
the dependence on the magnitude of the momentum. This leaves us with a function E{P) 
of a single variable. In this section we will describe the three experimental constraints which 
E{P) must satisfy if no new physics is introduced apart from the modification of the neutrino 
dispersion relation. 

4.1 Neutrino superluminaUty 

It is most likely that neutrinos are never superluminal. It is now clear that the superluminal- 
ity claimed by OPERA was due to experimental error [2]. However, in this demonstration of 
the use of neutrino splitting constraints we are attempting to understand general features of 
models which explain OPERA's neutrino superluminality claim. Therefore we will impose 
that neutrinos are superluminal in the region suggested by OPERA. This is certainly not 
a well-defined mission. It may well be that the neutrino velocity is a function of energy 
which is reasonably constant but has very deep subluminal troughs which are so thin that 
no OPERA neutrinos have hit them, and so they have not been observed. The existence of 
such troughs is very difficult to falsify, as they can be made as thin as one likes, and can in 
fact lead to a consistent dispersion relation, although one would need about 10 troughs in 
order to do this. Therefore we will assume that no such deep, narrow troughs in the neutrino 
velocity as a function of energy exist. For simplicity we will chose a constant velocity as a 
function of energy in the regime probed by OPERA's short extraction run. It will be clear 
that our result holds even if this constant is the minimum value allowed by the data, or even 
two or three standard deviations less. However we stress that our result is not qualitatively 
changed by any energy dependence of the dispersion relation, so long as the velocity function 
does not include at least 10 deep and narrow troughs. 

The group velocity of a neutrino is given by the usual formula from classical wave me- 



chanics 



dE{P) 
dP 
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and so this superluminal velocity corresponds to the slope of E{P), which is slightly greater 
than unity, throughout the regime probed by this second run. 



4.2 Cohen- Glashow Bremsstrahlung and ICARUS 

As explained above, if a neutrino requires more energy to arrive at a given momentum than 
two electrons, then phase space will be available for the neutrino to emit e~ pairs. Each 
such emission reduces the neutrino energy by a factor of order 1 [3J , and so if most neutrinos in 
any energy range experience such an emission, the spectrum will be distorted significantly. 
The strongest bound however comes from the high energy e~ pairs themselves. The 
linear dispersion relation considered by Cohen and Glashow would have resulted in millions 
of such pairs being detected by ICARUS, but in fact none were seen |29]. This fact was 
used by ICARUS to derive a bound on the neutrino superluminality. They assumed the 
linear dispersion relation of Cohen and Glashow, however this derivation may be applied 



relation to the function E{P) at other values of P. Therefore, the conservative bound that 
they give of 



for neutrinos near 20 GeV may straightforwardly be applied to the regime of neutrinos 
between 40-60 GeV. In fact, the real bound in this regime is stronger, due to the higher 
energies which far outweigh the fact that one is only using about 10 percent of the total 
neutrinos in deriving the bound. But this will be more than sufficient. 

Summarizing, we arrive at our second constraint. For neutrino momenta in the range of 
40-60 MeV, E < {I + IQ-^)P. 

4.3 Neutrino splitting 

The constraints from superluminality and from ICARUS may well be reconciled by a model 
in which neutrinos are subluminal at low energies and superluminal at high energies [TU]. 
In such models, E{P) stays below the diagonal line ii^ < P by at first having a slope 

■^Strictly speaking, this is only true for decays in which ah momentum is transfered into the electrons, in 
which case beta decay constraints on the neutrino dispersion relation yield a sufficiently low energy to justify 
our approximation. This can be extended to a sufficient fraction of phase space so long as superluminality 
stays below the 1 percent level. More generally, the phase space for this decay is dominated by values of 
P' in the OPERA energy range. In this regime the slope of E is greater than unity which means that the 
Cohen-Glashow process is kinematically allowed and likely sufficient alone to exclude these models. 



individually to neutrinos in any given momentum 
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of less than unity below OPERA energies, falling below the diagonal, and having a slope 
greater than unity at OPERA energies without ever climbing back to E = P. Before 
reaching E = P, above the energies probed by OPERA, the superluminality can be turned 
off with E < P at least up to 400 GeV, so as to be consistent with NOMAD. The region 
between 1 GeV and 10 GeV is constrained by MINOS [31], which does not convincingly 
demonstrate superluminal neutrinos but does convincingly demonstrate that neutrinos are 
not too subluminal. However, one can always compensate for this constraint by forcing 
neutrinos to be ever more subluminal at low energies. The only important constraint is 
that E{0) be smaller than 1 eV, so as to be consistent with bounds from beta decay. Thus 
a very concave E{P) may explain OPERA's observed neutrino velocity in a way which is 
compatible with ICARUS' lack of observed high energy electron positron pairs. 

The problem with this scenario is neutrino splitting [27], which as has been emphasized 
in Refs. [28] strongly constrains the function E{P). The constraints in those papers do 
not directly apply here, as we are allowing functions of p which can easily accommodate 
constraints from SN1987A, and we are not assuming any particular form of the dispersion 
relations. However, we will see that the applicable constraints are sufficient for our goals. 

In this process a neutrino with an initial momentum P and energy E{P) decays into 
2 neutrinos and 1 antineutrino with a total momentum P and total energy E{P). This is 
possible kinematically only if the function E{P) is concave or more precisely, as we saw in 
Sec. [21 if the group velocity exceeds the phase velocity 

in a large enough window of momenta beneath that of the incoming neutrino. 

Once the splitting process is allowed kinematically, it will slow down the neutrinos. The 
neutrinos lose an order one fraction of their momenta in each reaction. Like the Cohen- 
Glashow process, this will ruin the agreement between the neutrino energy spectra and that 
predicted by simulations in Ref. [12]. In particular it will lead to a maximum neutrino 
energy, but recall that we will only impose superluminality in the energy range in which 
OPERA has measured neutrino velocities in its short extraction run, and so we only have 
access to the dispersion relations in this regime. Therefore this is only in contradiction with 
OPERA if the energy spectrum is cut off inside of this energy regime, whereas a cutoff in 
the long high energy tail cannot be predicted. Finally, each of these processes will cause 
some deflections to the neutrino beam, making it less likely to hit the OPERA detector and 
therefore reducing the flux, again leading to tension with Ref. [T2]. In Sec. E]we will calculate 
the corresponding decay rate for OPERA neutrinos for dispersion relations consistent with 
constraints from neutrino velocity and ICARUS, and we will find that it is dramatically 
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inconsistent with the observed neutrino spectrum. 

5 opera's measurement of the neutrino velocity 

OPERA'S neutrinos are sent to Italy from the CNGS experiment at CERN. Most of them 
are the decay products of pions, which are created when 400 GeV protons ejected from the 
SPS in France hit a graphite target in Switzerland. The fact that the original protons have 
about 400 GeV each suggests that the neutrino energy distribution will have a tail extending 
to about 400 GeV. As the effects discussed in this note become more important at higher 
energy, for example with Cohen and Glashow predicting an energy loss rate proportional to 

[2j, it seems likely that future OPERA runs can lead to bounds which are much tighter 
than those in this note. This in particular applies to processes in which neutrinos emit /i /2 
pairs, 7r°'s and 7's. 

However during OPERA's 2009, 2010 and 2011 runs, due to the fact that it was using 
CNGS neutrino extractions which were 20 times longer than the effect for which it was 
searching, velocity measurements were statistical and indirect, and so could give rather little 
information about the energy dependence on the neutrino velocity. The results were most 
consistent with an energy independent distribution, but within two standard deviations could 
be fit by almost any distribution. 

The statistical error due to the long neutrino extraction time was essentially eliminated 
when OPERA received short pulses between October 22nd and November 6th of this year, 
allowing it to determine the velocities of individual neutrinos. However, given the short 
time dedicated to this experiment before SPS needed to begin accelerating lead ions, only 
20 acceptable events were produced. These events all yielded measured fractional super- 
luminalities between 1.7 and 3.9 x 10~^. This was interpreted as a more than 6 a signal 
for neutrino superluminality, but what does it tell us about the energy dependence of this 
superluminality? 

To determine this, first of all, one needs to consider the energy distribution of the neu- 
trinos seen by OPERA. This is shown in Fig. 10 of Ref. [T2], compared with the expected 
results from a simulation. As one can see, in all energy regimes the observed neutrino spec- 
trum is very close to that predicted, therefore effects like neutrino splitting and Cohen and 
Glashow's electron-positron pair production must be quite small throughout this range. In 
particular there is a long tail at high energies. While only 0.6% of the neutrinos from CNGS 
have energies about 100 GeV, the fact that higher energy neutrinos interact more strongly 
means that these account for 4% of all OPERA neutrinos [12] . 
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In what energy range do the 20 observed events he? It is impossible to know, since 14 
events occurred outside of the detector. However, based on the distribution in Fig. 10 of 
Ref. [12], one can determine statistically where they are likely to have occurred. As we are 
only interested in establishing a bound, we can safely use a very conservative estimate, and 
state that they all occurred between 10 and 50 GeV. What does the distribution of time delays 
then tell us about the energy dependence of the fractional superluminality e = {y — c)/c! 
Again, an energy independent value is favored, but a linear or reciprocal relation cannot be 
excluded. For simplicity we will use a constant distribution at the favored value 

6 = 2.5x10"^ (5.1) 

However we note that the splitting rate is proportional to e^, and so even if we used a 
constant value as disfavored as 10^^, this would only decrease the splitting rate by a factor 
of 20 which we will see is not sufficient to establish consistency. Moreover, a function which 
oscillates but rests within this regime, such as a linear or quadratic momentum-dependent 
fit, would lead to more splitting than a constant choice at the minimum value 10~^, therefore 
our analysis will also apply to these more general polynomial dependences. 

Now we are ready to assemble the superluminality constraint from Subsec. 14.11 with the 
ICARUS bound from Subsec. 1121 The former states that for 10 GeV < P < 50 GeV the 
slope is 1 + 2.5 x 10~^. The second states that during this entire interval it lies below the 
diagonal E' = P or if it surpasses the diagonal, it only surpasses by lO^^P. Of course this 
lO^^P is irrelevant, if indeed ^(50 GeV) = (1 + 10^^)50 GeV, then E(49.98 GeV) = 49.98 
GeV and so one may simply reduce the upper limit on P from 50 GeV to 49.98 GeV and 
conclude that E < P everywhere. For simplicity we will state 

E{P) <P In the range 10 GeV < P < 50 GeV. (5.2) 

So what happens for P < 10 GeV? It does not matter. We can consider neutrino splittings 
in which all three products have more than 10 GeV. The true rate, including all neutrinos, 
can only be larger, and so this establishes a lower bound on the decay rate. The decay rate of 
a 50 GeV neutrino to 3 neutrinos at energies between 10 GeV and 30 GeV can be determined 
using the dispersion relation described here. Notice that this dispersion relation depends on 
only two things. First, its slope is determined by superluminality which is measured directly 
by the OPERA experiment. Next, its height receives only an upper bound, the E = P 
line given by the constraint from ICARUS. Therefore in principle we have not determined 
the dispersion relation, it may always be shifted to lower energies. However such a shift 
will only make the dispersion relation more concave, and so increase the neutrino splitting 
decay rate. The most conservative calculation of the decay rate corresponds to the case 
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10 GeV 



50 GeV 



P (GeV) 



Figure 1: The dispersion relation. The velocity measurement at OPERA implies a slope of 
between 1 + 1.7 x 10"^ and 1 + 3.9 x 10"^ when 10 GeV< P <50 GeV. 

i?(50 GeV) = 50 GeV. In the next subsection we will determine the decay rate for this case 
and see that indeed neutrinos almost never arrive at OPERA from CNGS without splitting, 
in grave contradiction with Fig. 10 of Ref. |12] . 

6 Splitting decay rate 

We assume that in the range between 10 and 50 GeV the neutrinos can be described by 
the effective dispersion relation in Fig. [T] 



E{P) = {l + e){P-B) 



(6.1) 



where 



-5 



B 



e 



50 GeV ~ 1.2 



MeV. 



(6.2) 



2.5 X 10 



1 + e 



This is consistent with observations, indeed 



1. In the range between 10 and 50 GeV, the group velocity is 



dE 



(6.3) 



1 + 2.5 X 10 



2. At 50 GeV the neutrino 4-momentum is light-like {i.e. i?(50 GeV) = 50 GeV). 
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Finally, we have set the neutrino masses to zero, since they are negligible at these energies. 

The Lorentz violating kinetic term that can reproduce the effective dispersion relation 
dSH) is 

£fc = z/(po + B)^% + (1 + e)vpa'v (6.4) 
where e, B are the constants defined in (16.21) . 



6.1 Notation 

We are interested in the splitting process 

u{p) ^ u{p') + u{k) + u{k'). (6.5) 

Let us briefly clarify our notation: we use lowercase letters for the 4-momenta p, p', k and k' 
and uppercase letters for the magnitudes of the spatial vectors: P = |p| . The scalar product 
between two 4- vectors as p and k will be denoted by p-fc, while if only the spatial components 
of the momentum are involved we will write p-k. Ep will denote the zero component of the 
4- momentum p and so p ■ p = — Ipl"^ = — P^. 



6.2 Cross Section 

We will now calculate the decay rate for our neutral current process, neutrino splitting. The 
coupling between the neutrinos and the Z gauge boson is 



9 



2cosd 



w 



(6.6) 



The multiplicity factor of the relevant Feynman diagram is 2, and the corresponding ampli- 
tude is given by 

Y.\M\' = Y^ml \{u{p)^,u{p')) {u{k)rv{k'))\' . (6.7) 

spin spin 

In the Lorentz- violating theory in which the kinetic term is given by the (16. 4p . the complete- 
ness relations are 

u'{p)u\p) = v'{p)v'{p) = f (6.8) 

spin 

where 

p={p^ + B,{l + e)pi). (6.9) 
Similar notation has been introduced in Refs. j32[ 133] in this context. 
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After some straightforward algebraic manipulations, we find 

= 128G|(P-^')(P'-^) (6-10) 



spm 

and the decay width 



1 r dV f d-'k r d-'k' l,.,,o, ,4,4, , , ,n 
(2TTfEj J Et J 2Ei,,"^ I w I- I (2j)5£)^ 

(6.11) 

We will restrict the domain of integration to energies between 10 and 50 GeV, since 
OPERA only indicates that fl6.ip holds in this range. This means that the value F that we 
will obtain is only a lower bound; however we will see that it is sufficient for our purposes. 
More or less following the phase space integration strategy of Ref. [9] we find 



dV f(i k 



Epi J Ek 



(Tkiv ■ k'W ■ k) 5{k' ■ k') 6\p-p' -k-k') (6.12) 

E,j>10 



^ ^ ' (j) . (j) — p' — k))(p' ■ k) 6{\p ~ p' — fcp) 



Epi I Eh 



Ep-E^,-Ek>10 



where |p — — A;^ = {p — p' ~ k) ■ {p — p' — k). 

It is convenient now to define 9i as the angle between k and p—p' and perform the change 
of variables x = cos6'i. Similarly, we let 62 be the angle between p and p' and y = cos02- 
We will need the relations 



J, _ pj. (P - P'y)^ + p'Vi^'Vi^' cos ^g^^g^ 

\p-p'\ 

and 

p' ■ k = p ■ k — {p — p') ■ k = p ■ k — \p — p'\Kx. (6.14) 

We can write the last delta function as 6{f{x)) where the function f(x) has a single root 

\p - pf -{P-P' + 2Bf + 2k{P -P' + 2B) 

xo = (6.15) 

2K\p-p'\ 

and its first derivative reads 

f'{x) = 2{l + efk\p-p'\. (6.16) 
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After integrating over x, Eq. f l6.12p reads 

d^p' f KdKdipi {p- {p — p' — k)){p' ■ k) 



Epi J Ek 2\p — p'\ 



Ep~Ep,-Ek>W,x=xo 



The domain of integration is determined by the conditions — 1 < Xq < 1- Moreover, our 
phase space is hmited to the region in which all the neutrinos have energy between 10 and 
50 GeV, yielding the additional conditions 

P-P'-K>WGeV (6.17) 
K >10 GeV. (6.18) 

Notice that in this case we have neglected the contributions to the energy proportional to 
e. Indeed, following the same argument used in the previous section, one can easily see that 
they are equivalent to an irrelevant shift of the region of validity of our effective dispersion 
relation. 

From the conditions on Xq we obtain 

P-P' + 2B-\p-p'\ ^ P-P' + 2B+\p-p'\ 

^ <K< ^ (6.19) 

2BiB + P-P') , , 

^ > 1 - ^ p.p -■ (6-20) 

If the condition fl6.20p is not satisfied, Xq < 1 implies 

which is in contradiction with the condition f l6.17p . Since it easy to see that fl6.17p and (I6.18P 
are stronger bounds than (16.191) our domain of integration is given by 

10 GeV < K < P - P' - 10 GeV (6.22) 

2BiB + P-P') 
y > 1 ^ • (6.23) 

The conditions IK22^ also imply that 10 GeV < P-P'- 10 GeV ^ P' < P-20 GeV, which 
yields the last bound. Choosing P = 50 GeV and performing the integration we obtain the 
decay rate 

r = 3.8 X 10"^^ GeV (6.24) 

and the decay length 

he 

£ = — = 0.52 km. (6.25) 
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This is much shorter than the observed 730 km distance that neutrinos travel from CNGS 
to OPERA. Therefore such models are strongly excluded, because the observed neutrino flux 
would not fit simulations as well as seen in Fig. 10 of Ref. [12]. As in the case of Ref. |3] the 
decay rate at a fixed energy is proportional to the cube of the difference between the electron 
and neutrino dispersion relations, therefore a suppression of this process requires that these 
two energies be at least an order of magnitude closer. As the slope at most momenta P is 
fixed directly by the OPERA velocity measurement, this requires that the neutrino energy 
curve in Fig. [T] by modified by including narrow corrections in which the neutrino energy 
rapidly approaches the electron energy, corresponding to a massively subluminal neutrino in a 
very narrow momentum band. To reduce the difference between the two dispersion relations 
by an order of magnitude, one requires at least 10 such narrow bands. In a polynomial fit 
to the dispersion relation, this requires a polynomial of order much greater than 10. 

7 Conclusion 

In this note we have seen that the neutrino splitting process u ^ 2u + u places strong 
restrictions on the form of the neutrino dispersion relation E{P). The result is effectively 
that E{P) will be linear or convex, or more precisely that the group velocity will never 
exceed the phase velocity. This is a quantum analogue of the modulational instability in 
classical wavemechanicqj in which nonlinearities (corresponding to the self-interactions in 
the quantum theory) cause a plane wave to decompose into pulses. As in the quantum case, 
this occurs when nonlinearities in the dispersion relation yield a group velocity greater than 
the phase velocity in a certain range of the space of momenta. However in the classical case 
the number of excitations does not change, the instability serves only to drive the Fourier 
components of the waveform away from the original planewave momentum. 

7.1 Applications for neutrino splitting 

As the splitting rate depends only upon the neutrino sector, it may be used to place robust 
constraints on a theory, independently of the dispersion relations of the other particles. It 
therefore has numerous applications. In this note we have seen that it implies that a change 
to the neutrino dispersion relation alone, even with an arbitrary density dependence, is not 
sufficient for consistency with OPERA and ICARUS results. 

The concavity bounds may also be used to place upper bounds on superluminality at 

"^For a review of various manifestations of these instabilities in classical physics, see for example Ref. 134' . 
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high energies, for example at ICECUBE. Indeed, using the fact that 400 TeV neutrinos 
are depleted by at most a factor of 2 with respect to simulations [E], one may reverse the 
argument of Ref. to exclude the e ^ E and e ~ E'^ behaviors which historically have 
been expected from effective field theories. 

7.2 Experimental signatures 

In closing, we would like to mention that if indeed the neutrino dispersion relation does 
depend on density, but that this is consistent with the above constraints because for example 
charged lepton dispersion relations have a similar dependence, then there are several kinds 
of experimental signatures for which one may search. 

Besides the neutrino, there are two particles that may travel long distances, of order 
a kilometer, through solid rock. These are the and the 7r°. Millions of them, created 
by interactions with cosmic rays, have been detected in deep detectors such as MINOS. 
If their dispersion relations are modified with respect to those of other particles, one may 
look for signatures of this change in their decay rates inside of various media or even direct 
measurements of their travel time from the surface to underground detectors. 

However it may be that the dispersion relations of all particles have the same density 
dependence, then what signature is left? Recall that in 1979 FermiLab [35] was able to place 
strong bounds on neutrino superluminality, of order 5 x 10^^, by racing neutrinos and muons 
a kilometer through rock. Significant error was introduced due to uncertainties involving 
the time lost by a muon traveling through rock. However this experiment was intrinsically 
more precise than experiments like OPERA and MINOS because there was no reliance on 
an external clock, one merely measured the time difference the arrivals of the muons and 
neutrinos at the same location. This is why a long baseline was not necessary, one kilometer 
already gave a precision comparable to MINOS' 734 km experiment 25 years later. 

One can do a very similar experiment racing neutrinos passing through rock and neutrinos 
beside them passing through a tunnel. If one is lucky enough to receive neutrinos that 
passed both through the rock and through the tunnel in the same extraction, then one can 
compare their arrival times and thus measure their relative velocities very accurately. Such 
conincidences will be much more frequent at MINOS' near detector after the energy increase 
of its neutrino beam scheduled to take place in two years' time. 
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